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Abstract 

In inclusive P-wave charmonium decay, to cancel infrared divergence in the 
colour singlet contribution requires the inclusion of the colour octet which becomes 
degenerate with the former in the infrared limit. In the corresponding exclusive de- 
cay, such an infrared divergence does not exist. On this ground, it becomes doubtful 
whether the colour octet is needed in exclusive reactions. A contradiction in the un- 
derlying picture, however, would result once all strong decay channels are summed. 
We provide an answer to this question in support of the colour octet with theoretical 
arguments as well as an explicit calculation of xj — ► PP- 
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1 Introduction 



Within the standard hard scattering approach (SHSA) of Brodsky and Lepage U for 
exclusive reactions and the later modified version (MHSA) of Botts, Li and Sterman ||], |3j , 
higher Fock states are usually suppressed due to the necessary hard momentum exchange 
for ensuring that the constituents of each outgoing hadrons are almost collinear. As a 
consequence, only the knowledge of valence Fock states is necessary. However, there are 
exceptions. For example, this might not be true when there are suppression mechanisms 
at work due to helicity conservation, flavor symmetry, higher- wave wavefunctions etc.. 
Indeed, it has been shown in inclusive P-wave charmonium decay [|], [5| that the next 
higher Fock state, the so-called colour octet, where the constituents of the charmonium 
consist of the usual charm-anticharm pair and an additional constituent gluon, is needed 
not only for an infrared safe inclusive decay rate but also because they are of the same 
order within a non-relativistic QCD velocity expansion formalism. 

Turning to exclusive processes, the same octet mechanism should also be at work in xj 
decays. However without a manifestation of an infrared divergence, the necessity of the 
colour octet was never there. In this paper, we shall try to clarify this issue. The paper is 
organized as follows. In the next section, we will describe the situation in inclusive decay 
and provide the main motivation for what are to follow. In Sec. [3], the role of the orbital 
angular momentum of the quarkonium and its effect on the annihilation are explained in 
a language consistent with the calculation scheme we use below. Then power counting 
arguments for the colour octet is provided in Sec. |j. The remaining sections are dedicated 
to explicit calculations. First we will show in Sec. |^ that the colour singlet contribution 
in xj — > PP within the MHSA is not sufficient, and then the method of calculating the 
octet contribution within the SHSA will be explained in Sec. || The results of combining 
the singlet and octet contribution within the latter framework will be given and compared 
with experiments. 

2 Is Colour Octet Necessary In Exclusive Decay? 

In the early 70 's, inclusive decay of charmonium was studied by Barbieri et al || [7j 
in a series of papers. The family of quarkonia, which consists of members made up of 
bound state of heavy quark-antiquark pair, are special because they can be treated non- 
relativistically. These systems decay strongly via annihilation into gluons. The P-wave 
charmonia, which are our central interest in this paper, decay through two gluons at 
leading order in a s (Fig. [I] (a)). In the rest frame, two gluons depart back-to-back with 
their kinematics completely fixed. The decay probability amplitudes are therefore regular, 
free from any divergence as a result. This is the case for xo and X2 at leading order. For xi, 
the same decay into two massless spin-1 gluons is forbidden. The leading process in this 
case is the decay into a gluon and a quark-antiquark pair. A three-body final state is less 
restrictive and has more kinematic freedom. In particular, the quark and antiquark can 
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Figure 1: The leading decay process for (a) xo and Xi an d (b) Xi- A possible decay of 
the higher Fock state, the colour octet, of the xj is shown in (c). In the infrared limit of 
the outgoing gluon, (b) and (c) are degenerate. 

now depart near back-to-back leaving the gluon to take the remaining energy-momentum 
and thus ensuring the latter two are conserved. This means that there is no lower limit 
to how soft the gluon can be. It is well known that there is an infrared divergence in the 
tree graph as shown in Fig. [1] (b), that is, provided that all external legs are put on the 
mass-shell. In view of the divergence, it was decided in || [7] that since the cc was really 
coming from a bound state so it was perhaps more appropriate for the heavy quarks to 
be off-shell. This off-shellness must correspond to the bound state energy e. With this as 
a cutoff, the inclusive decay width acquires a logarithmic e dependence 

The same problem actually occurs also in the inclusive decay of Xo an d X2 but at next- 
to-leading order. This form of the logarithmic dependence on the binding energy of the 
width was left alone until the 90 's. The emergence of the colour octet higher Fock state of 
charmonium and the realization of its importance to P-wave charmonium enable finally 
the cancellation of this infrared divergence @, |5j. In the infrared limit of the emitted gluon, 
the singlet decay graph becomes degenerate with the annihilation of a cc in a colour octet 
state into qq with an accompanying soft spectator gluon. This is illustrated in Fig. |1] 
(c). The infrared divergence cancels between Fig. [I] (b) and (c). Thus in inclusive decay 
because of infrared divergence, it is incorrect to separate the valence colour singlet state 
from the higher colour octet state for a P-wave charmonium. Both have to be included to 
ensure infrared safe calculations. The same, however, cannot be said in exclusive decay. 
By stipulating what final states should be observed, the infrared problem is removed. 
The kinematics of decaying from a P-wave charmonium bound state into other bound 
states does not allow any infrared divergence. Since no problem was encountered, it was 
traditional to calculate the partial decay widths using only the colour singlet valence 



Fock state |8], [], [10], [□], [L2| and claims were made that these calculations could give 
results in agreement with experimental measurements. In view of the development in 
inclusive decay, it is unavoidable to ask the question whether colour octet is needed if 
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one concentrates on just a few decay channels. If the answer is no to each and every 
channel, then there would be a contradiction since the partial width of each individual 
channel could in principle be worked out separately and the sum of them must agree with 
the total hadronic width. In practice, if the requirement of the colour octet does not 
reveal itself in at least one of the few-body decay channels, it would be hard to find it 
in the channels with many final products. To confirm this contradiction would then be 
very daunting. However it is hard to imagine how the colour octet would be important in 
many-body decays but not in the few-body channels since the octet comes in already at 
the leading order of xi an d at the next-to-leading order of xo and X2 in inclusive process. 
With the lack of the same infrared divergence, one needs a more compelling reason than 
the less tangible claim of a contradiction to introduce the colour octet in exclusive decay. 
We will provide an argument to this effect in the following two sections. 

3 Angular Momentum Barrier To Annihilation 

Charmonium decays strongly via annihilation into gluons. The heavy mass of the cc 
system means that the heavy quark and antiquark tend to come close to each other when 
the annihilation is taking place. The probability of small spatial separation of the cc will 
then affects the likelihood of the annihilation to take place. For a S-wave charmonium, 
the probability of spatial distribution is largest at small distance so it is favorable for the 
annihilation. But for a P-wave, the L = 1 orbital angular momentum of the cc forces 
them apart and the probability for the heavy fermions to be near each other tends to be 
very small and, in fact, vanishes at zero separation. For this reason, the annihilation of a 
P-wave charmonium into gluons is suppressed by this angular momentum barrier. 

To get an idea of the actual size of the suppression, one only has to examine the 
charmonium wavefunction, which contains all the information of the spatial distribution 
of the bound state system. In the case of a S-wave, annihilation occurs at a length 
scale / ~ 1/M, the inverse of the charmonium mass. The probability amplitude for the 
annihilation then depends on ips{l ~ — V's^ ~ 0) the wavefunction at the origin. 

As for a P-wave, because of angular momentum ipp(0) = 0, so the product of I and the 
first derivative of the wavefunction has to be used instead ipp(l ~ 1/M) ~ I if)'p(l ~ 0). 
Fourier transforming these wavefunctions into momentum space, we have 

S-wave: ?As(0) — > ^ps(k) 

P-wave: I r(/ P (0) — ^p(^) , 

where k is the internal transverse momenta of the quarkonium of the order of a few 
hundreds MeV. It is now obvious that the suppression in the P-wave wavefunction is of 
the form 1/M which is substantial for a heavy charmonium to decay through annihilation 
into gluons. 
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4 Power Counting in the Charmonium Mass in The 
Decay Probability Amplitude 

In the last section, we showed that P-wave charmonium annihilation was unfavorable in 
comparison to S-wave because of the angular momentum barrier. To answer the question 
we posed earlier, it remains to show how this fact brings about the need of the colour 
octet. To do that we need to examine the large scale dependence of the decay probability 
amplitude or more specifically to perform power counting using that scale. The simplest 
scheme to work out a decay amplitude is via the hard scattering scheme already mentioned 
at the beginning. In this SHSA method, the decay amplitudes for a charmonium C = 
J/i[),Xj, ■ ■ ■ to decay into a pair of pseudo-scalars V = 7T,K,... etc. and a proton- 
antiproton pair pp have the form 



Mc—>w ~ fc4>c ® fv<Pv ® fv4>v ® T H (2) 



and 



Mc—,pp~fc4 > c < 3fp4>p < 8fp4>p®T H , (3) 

respectively. In Eqs. (Q) and (0), the amplitudes are each given by a convolution over 
probability amplitudes 0's and the hard perturbative part Th ■ The partial decay width in 
each case is given by V ~ \Ai\ 2 /M. So M. is of mass dimension one. In these amplitudes, 
only the decay constants /'s and T# carry a mass dimension of some power. These must 
altogether make up for the right mass dimension of M.. Each product f<j) is nothing but 
the internal transverse momenta integrated light-cone wavefunction 

f(f)(x) = J J[ -^-^ ip{x;k ±1 ,. . . ,k ±N ^) (4) 

where Q is some large cutoff scale, given approximately by the large scale of the hard 
process in question, for the N — 1 internal momentum integrals over the N particle 
wavefunction ip, and x represents collectively (x±, . . . ,x^) the light-cone fractions. The 
normalization of the wavefunction is given by 

II ^ Xi ~(2^ \ip(x;k ±1 ,...,k ±N _ 1 )\ 2 = P . (5) 

Here P is the probability of this particular hadron to be in this state with iV constituents 
and is of course a dimensionless number. Using Eqs. @) and (|^), the mass dimension of 
the decay constant of any hadrons can be deduced. In this paper, we restrict our consid- 
erations to only charmonium states with up to three constituents for any charmonium. 
In Table [I], the mass dimensions of the decay constants together with the number of con- 
stituents of several hadrons are listed. The odd ones in the table are the decay constants 
of the colour singlet Xj\ which has a dimension of two although it is coming from a 
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2-particle wavefunction, and the colour octet J/ip^ 8 \ which is a 3-particle wavefunction 
but has a dimension of three. The reason for this is that as we have already discussed, 
the P-wave wavefunction vanishes at the origin so the product of the first derivative of 
the wavefunction at the origin and the annihilation size I must be used instead. This 
has the effect of increasing the mass dimension of the decay constant by one. In the case 
of J/ip^ 8 \ the fact that J/tf> is a 1 charmonium excludes the possibility of the leading 
colour octet contribution to be in a S-wave three-particle wavefunction, which must be 
found in a three-particle P-wave state instead. Note that for the three-constituent colour 
octet state of the J /if), the cc pair itself can be in 3 Pj or 1 So as is frequently considered in 
inclusive J /if) productions jl3[ [14], 15], HU, but the full wavefunction must have an orbital 
angular momentum of one. In Table |2], the orbital angular momentum L% and spin S C c 
of the cc pair, the orbital angular momentum between the pair and the constituent gluon 
L c g c , the number of constituent gluon N g , the parity V and C-parity C of each state of the 
J /if) and xj u P to three constituents are tabulated. This shows that the 1 S colour octet 
state of the J/ if), although has no orbital angular momentum in the cc pair, must have 
one unit of this momentum between the pair and the constituent gluon otherwise this 
state cannot be in the required 1 state of the J /if). So although the heavy quark pair 
is in an S-wave, the three-particle wavefunction is effectively in a P-wave. That is why 
the decay constants of the two three-constituent colour octet states of the J /if), usually 
labelled as 3 Pj and 1 S'o, are both of mass dimension three. Therefore the mass dimension 
indicated in Table [I] for J/ipw is true for both states and L = L\ + L^ c in these two cases. 
In Table |2], two possibilities of the potentially possible three-constituent octet states of 
the J /if) with the heavy quark pair in 3 Si, but which in reality do not exist hence the 
attached "?" marks, are shown. The purpose of these two entries in the table is to show 
how using the remaining available degree of freedom L c g = or 1, it is still not possible 
to change the states into a 1 because L g c only affects V and not C. In the case of the 
octet 1 S'o, it is the L c g c = 1 which permits this possibility to exist. 

Since the decay process of C has only one large scale, that is the charmonium mass 
M, Th must contain sufficient power of this scale to make up for the mass dimension of 
M.. So in the following decay processes, the power counting in or the power dependence 
on Mj/^ ~ M XJ ~ M of the charmonium in the colour singlet M.^ and octet Ai^ 
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amplitude with up to three constituents in the respective charmonia goes as follows Q 



As) , , 



'X.J 1V± XJ 

AS) , 

"2— ~ ^#(^) 2 ~^ (12) 



M 8 ^- - ^4h^) 2 ~4- (13) 



M2 / M 5 

Eqs. (||), (P), flT2"D and (|TB| ) show that the probability amplitude of the colour singlet and 
octet of x.j decay depend on the same power of M. The contribution from the colour octet 
higher Fock state of the P-wave charmonium is not suppressed at all as naively expected 
in comparison with the singlet contribution. However, this proves to be true in the case 
of the S-wave J /if). Eqs. fl6|), ©, (0) an d (HH) show that the colour octet contribution 
to the amplitude of the J/if> decay is indeed suppressed by 1/M 2 (this suppression would 
have been only 1/M if the leading three-particle wavefunctions of the colour octet states 
of J/ if) were not forced to be in a P-wave when all three constituents have been taken 
into account, as such there is an additional factor of 1/M). The power of M counting 
of the octet J/ if) decay amplitudes in Eqs. (|7|) and ( |TTD are applicable for the cc in both 
the 3 Pj and 1 5'o state. Comparing the colour singlet decays of the J /if) in Eqs. (H) and 

3 Note that our power counting is done by taking the charmonium mass M to be the largest and most 
dominant scale in the process and the ratios of any quantities with any others are much less important 
than those of M with these quantities. In this sense, one can forget about the actual magnitude of the 
other quantities with a mass dimension, such as the decay constants, with the understanding that they 
are small in front of M. This permits us to concentrate only on the power of M in the amplitudes above. 

4 We have, for the sake of simplicity, ignored the dependence of the charmonium decay constants 
on the large scale M because this is at most logarithmic and is therefore much weaker than the power 
dependence. 
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( fLOf) with the \ J decays in Eqs. (||) and (|i~2|) shows that the J / ip amplitude is larger by a 
power of M in each case. The reason for this is due precisely to the fact that the P-wave 
charmonium wavefunction is suppressed by 1/M because of the L = 1 angular momentum 
as explained in the previous section. Without this angular momentum suppression, the 
colour octet contribution in \j decays could be neglected. But as it is, provided that 
there is no other form of significant suppression, which is indeed true fT7| , |TB| , [T9|] , colour 
octet contribution must be included. 



5 Colour Singlet Contribution in the Modified Hard 
Scattering Approach 

The previous sections provided the arguments for the inclusion of the colour octet in xj 
decay. The decay into pseudo-scalars with the octet contribution has been considered 
in [|Tj], |I8"|| . In this section, we will show by explicit calculation, that the size of the 
colour singlet contribution to the decay into proton-antiproton is small in comparison 
with the experimental partial width. The scheme that we will use is the MHSA because 
of two advantages it has over the SHSA. The first is the renormalization scales of the a s 's 
are determined dynamically by the virtualities that flow in the diagrams that represent 
the decay process. The second is the problematic endpoint regions of the distribution 
amplitudes where the virtualities tend to approach Aqcd are suppressed by the inclusion 
of the Sudakov factor exp{— S}. The decay probability amplitude takes the form 

Mxj^pp ~ ^xj ® i>p ® tpp ® T H (®a a ) <g> exp{-S} . (14) 

The a s in T# is expressed in such a way to signify that the coupling, although hidden 
within the hard part, is part of the convolution and is not taken outside as a prefactor as in 
the SHSA. As can be seen, a proton wavefunction has to be used. This is potentially prob- 
lematic because of all the uncertainties surrounding the proton distribution amplitude. 
Although proton model wavefunctions exist via QCD sum rule derivation (see any of the 
refs. !| ^0], ^1|, [23|, [23j), none of them could describe the magnetic form factor below 



Q 2 ~ 50 GeV 2 in a self-consistent manner |25|, 26 1. Fortunately by relinquishing a formal 



derivation, one can turn to a phenomenological approach. Since the high momentum tail 
of the wavefunction is clearly not dominant in present day experimental accessible regions 
- Q 2 even up to 70-100 GeV 2 are still nowhere near the asymptotic region as far as the 
proton is concerned — the soft part of the wavefunction must still be dominant up to 
this range. Then by fitting the soft overlap of the initial and final proton wavefunction, 
or the so-called soft Feynman contribution to the magnetic form factor, to experimental 
constraints such as the J/ip — > VP decay width, proton valence quark distribution and 
the magnetic form factor measurements, a model proton wavefunction can be constructed 
|27| . This, by construction, works automatically in the Q 2 ~ 10 GeV 2 region, which is 



precisely the region of interest in this paper. We stress that the problem discussed here 
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concerning the proton wavefunction is strictly limited to this sector of the convolution in 
Eqs. (0) and (||). It does not affect the charmonium or the perturbative hard part because 
as in the problem of the proton electromagnetic form factor, it involves only the proton 
wavefunction. Once we used the phenomenologically constructed wavefunction and not 
those from, for example, sum rules in this sector, the problem is taken care of. 
The proton wavefunction is given by 

\P, + ) =^/[d*][d 2 fcJ {*? 2 3|<^)+^ 13 |^<^+> 

-(*?32 + *5si) |« 2 ^>} (15) 

where the decay constant f p and the distribution amplitude (pi 23 in the scalar function 

1 



123 



[x,k±) = V(x 1 ,x 2 ,x s ;k±i,k±2,k±- i 



8V3! 



/pOM^fo Vf) • ( 16 ) 



depend on the factorization scale /xp, and the two measures in square brackets are the 
usual ones 



3 3 

[da;] = n dxi 5(1 - x i 

8=1 1=1 



[d 2 fcj 



r 3^2 



3 3 

fld 2 k ±i 6®(£k ±i ). (17) 



The internal transverse momentum dependent part can best be expressed as a Gaussian 



Qp(x, fcj_) = (167T 



2\2 



XlX 2 X 3 



-cxp 



X j 



i=l 



;is) 



where a p = 0.75 GeV 1 is the transverse size parameter and f p (no) — 6.64 x 10 3 GeV 2 
at the reference scale /iq = 1.0 GeV. These values are obtained in the construction in 



27J ]. The proton distribution amplitude in terms of the Appell polynomials </>i 2 3's, the 



scale dependent expansion coefficients B^(^f) and the asymptotic distribution amplitude 

0as(^) is 



0?23 (X, Hf) 
0123(^^0) 



0As(^) 
(f>As{x) 



n=l 



(19) 

60xia;2a;3[l + 3a;i] . (20) 



Eq. 



shows the rather simple form of 0i 2 3 at the scale //q- For the scale dependent 



expression of B^(hf) and / p (//f) and those of the more general baryon wavefunctions, we 



refer the readers to our other paper |19| . 

Equipped with this proton wavefunction, we are ready to work out the perturbative 
hard part Th in MHSA. Because of C-parity conservation | 28fl , the colour singlet xj 
component can only annihilate into proton-antiproton through two gluons [I9|. The 
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Figure 2: Feynman graphs for the colour singlet \J decay into proton-antiproton. 

Feynman graphs are depicted in Fig. |2]. The three light quark lines in each case can be 
permuted further to give more graphs but these can be included by a symmetry factor. 
Basing on Fig. [| (a), T^(x, y, b, b') in transverse separation space b and b' of the proton- 
antiproton and in terms of 



C({x,y) = y 3 (l-2 Xl ) 2 - J 
C J 2 {x,y) = (-1) J_1 y 3 xi(x! - Vl ) 



(21) 
(22) 



worked out to be |19| 

fJ{x,y,b,b') = 



2 9 x/2 



9V3 



a a (tiK(t 2 K(t 3 ) tf2 ( 6 i ~ b i ~ b 3 + b' 3 



Xy^ 1) (V^(2m c )|6 3 |) ^ 1 \^l-x 1 )y 3 (2m c )\b 1 - 6&|) 



H^ 1 \^m(2m c )\b' 1 \) 
(x x +y + 1)(1 -Xi -y x ) 

^ 1) (^/(l-x 1 )(l-y 1 )(2m c )|6i|) 



C 2 J (x,i/) 



1 -xi -yi){2-x 1 -yt) \ 2-x x -y^ 



4 ^(^(l - y x ) + yi (l - ^Qmjjj) 
7T (xi + yi)(2 - -2/1) 

(cfov) + , + 2C Jt' v) — 

v {x 1 +y 1 ){2-x 1 -y 1 ) / 

2m c \b[\ Ci(x,y) K^Ux^l - Vl ) + Vl (l - xJmM 



yxiO- - Vi) + - JPi) (2 - x 1 - yi)(x! + y x ) J 
+ ((x, b) < ► (y, 6')) • (23) 

The f/o is the Hankel function, ITq and are the modified Bessel functions of the 
second kind. In the above and throughout this paper, we use the one-loop expression for 
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a s and set Aqcd = 220 MeV and m c = 1.5 GeV. The renormalization scales t\,t2 and £3 
in the product of three a s 's here are determined by the largest of either the virtualities 
in the neighborhood of the vertices associated with each a s or the smallest transverse 
separation b or b' of the constituents of the proton or antiproton in accordance to the 



method described in Q, F3. 26 



Convoluting this with the wavefunctions and the Sudakov factor gives the decay form 
factor 

_ i vmm^ f ldxUv] fj (x v b bf) 



#? 23 (x, b)V p 123 (y, b') + * p m {x, b)* p 321 (y, b') 
+ (^ 23 (x, b) + ¥ m ( x , 6))(*? 23 (y, 6') + fc&ufo, b>) 
+ (2^3)} ex V [-S(x,y,b,b',2m c )} (24) 

where oj = l/\/2, 1 for J = 1, 2 respectively, and the first derivative of the P-wave radial 
wavefunction of xj is well known |-Rp(0) | = 0.22 GeV 5 / 2 . The form of the Sudakov function 
S can be found in, for example, ref. [fl8|| . The Sudakov exponential factor performs a very 
important function of suppressing the regions where any of the virtualities U approaching 
the dangerous scale of Aq CD . It renders the calculation self-consistent in the sense that 
the size of the contribution from the soft region can never be the bulk of the hard part. 
With this self-consistency in place, the decay form factor is related to the decay width by 

r(xi - Pp) = p ^ Mp,M ^ - Y] \m\ a a 2 = p ^ Mp ' M ^ < \B P \\ (25) 

'' \QttM 3 ' A1A2A q^/i# 



XI " A's ""^"Xl 



3vrM v 



and 



r / v tmi _ gPjWy 1 y|u2 2_ Pp , s .(M p /M X2 )mg 2 



16vrM X2 5 Vg 



10ttM X2 



Since in the calculation proton is taken to be massless, phase space has been corrected by 
p p . s .(z) = VT^I^. 

Our results are shown in Table § together with the experimental measurements from 



Particle Data Group (PDG) |29| and the more recent BES Collaboration | 30| . The errors 
associated with the PDG results are shown but those associated with BES are not because 
they are too large. One can consult them in ref. |30| . The results reveal a large discrepancy 



between the colour singlet contributions and the experimental widths. Before drawing any 
definite conclusions from these results, we first show that our calculation within MHSA 
is perturbatively consistent. To do that, an artificial cutoff a c s rit ' in the value of the 
coupling was introduced so that the integrations in Eq. (Elf) include only regions where 
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Figure 3: Percentage of the colour singlet contribution to the xj — > PP decay width 
that comes from region of values of a s below some critical cutoff value af lt \ The solid 
and dot-dashed line are for xi an d X2, respectively. This shows that the bulk of the 
contributions in fact comes from a s < 0.6 region. 



Y[i=i a s(ti) < (a" 1 ') 3 was true. The value of a™*' was varied to determine how much of 
the final results came from hard regions. The percentage of the final widths are plotted 
against a™*' in Fig. |3| for both J = 1 and J = 2 charmonium. It shows clearly that by 
far the main contributions come from hard regions in which a s remains below 0.6. So 
the Sudakov factor performed as expected to keep the soft non-perturbative region under 
control and suppressed. 

To assess the uncertainties in the results, we have used the freedom from the uncer- 
tainty of the relation between R' p (0) and m c from fit to charmonium parameters [31) and 
from quarkonium potential model [^] and also the value of Aqcd to maximize our results 
to see if they can be raised up to the experimental results. We found that they maximized 
with values of R' p (0) = 0.194 GeV 5 / 2 and m c = 1.35 GeV and a larger A QCD = 0.25 GeV. 
The colour singlet contribution to the decay width remains for J = 1, 5-9 times and 
for J = 2, 2-4 times below the measured results. At this point, one may wonder that 
the uncertainties of the proton distribution amplitude or wavefunction may make up for 
the remaining difference. As we have already discussed above, many proton wavefunc- 
tions derived from QCD sum rules do not work at our present scale of interest. For this 
very reason, the particular one we used here was constructed phenomenologically. Thus 
the apparent many choices of wavefunction and the fact that using each and everyone of 
them may lead to very different results are illusionary. One simply cannot use these other 
wavefunctions for xj decays. As to the nucleon decay constant, other values do exist from 
QCD sum rules or from lattice. For example, the values from QCD sum rules is about 



5.0 x 10 3 GeV 2 fl2|, [2T| or from lattice QCD there exist the two values 2.9 or 6.6 x 10" 



GeV 2 0. Note that these are all smaller than ours, so using them will give even smaller 
results for the singlet contribution and not larger. The large experimental discrepancy of 
the data between PDG and BES cannot be explained but assuming that the true widths 
are somewhere in between or even near the BES results, it is clear that the colour singlet 
contribution alone cannot account for the experimental partial decay widths in agreement 
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with our theoretical arguments presented in Sec. [| and f|. 



6 Combining Colour Singlet and Octet Contribution 
in the Standard Hard Scattering Approach 

In the last section, we have shown that colour singlet contribution is not enough to account 
for the experimental measurements contradicting all earlier similar calculations done be- 
fore the emergence of the colour octet state in P-wave charmonium and the associated 
improved understanding of the role of the higher Fock state of the quarkonium. It re- 
mains to show that the colour singlet and octet combination do indeed provide reasonable 
agreement to experimental data. We will outline the method of dealing with the colour 
octet and describe our calculations. The combined results will be shown but the details 
are given elsewhere [TSI . 



The colour singlet calculation above was performed within the MHSA because of the 
two already mentioned advantages. The drawback of this scheme is, however, that it is 
more complicated to do a calculation. One has to include full wavefunctions and not only 
the distribution amplitudes. This means internal transverse momenta have to be kept 
in the hard perturbative part T#. For the colour singlet contribution to xj decay, this 
proves to be manageable. We have seen the transverse separation space version of Tjj in 
Eq. (^). It contains several hyper-geometric functions and is therefore non-trivial. In 
the case of the colour octet, we will see presently that the calculation involves many more 
diagrams. In addition, proton and antiproton each carries three valence quarks. All these 
add up to the dimensions of the numerical integrations and complicate the calculation 
much more. The previous calculation was done within the MHSA to avoid the ambiguity 
of having to choose a renormalization scale by hand for a s since the decay widths depend 
on the sixth power of the coupling a Q s which leaves too much freedom or there is a lack of 
constraint to arrive at results that agree with data. In view of the complexity of MHSA, 
we will keep the calculation simple and avoid tangling ourselves in MHSA in the high 
dimensional computations by retreating back to the SHSA. This will, of course, bring 
back all the shortcomings of the SHSA but in spite of these, the SHSA can often be used 
to make a reasonable first estimate of the results. 

As is now well known, the colour octet state of a P-wave charmonium includes, in 
addition to the usual cc, a constituent gluon. The presence of this in the initial state 
requires it to be properly "put away" in the final hadrons in order not to have a net 
colour. This colour neutralization can be achieved in two different ways. The first is to 
allow this gluon to enter directly into or be a constituent in the final wavefunction of 
one of the outgoing hadrons. An example is shown in Fig. |] (a). SHSA requires that no 
participant in the hard interaction to be disconnected from the rest as a consequence of 
the collinear approximation in SHSA. That is all constituents of any hadrons, be they part 
of an incoming or outgoing hadron, must be almost collinear. Therefore the connectedness 
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— \smsL 



(a) (b) 

Figure 4: Example of two ways to achieve colour neutralization: (a) the constituent gluon 
enters the nucleon as part of the constituent in the nucleon after exchanging a gluon with 
the other constituents, (b) this gluon ends in the perturbative hard part T#. 




Figure 5: Example of cc annihilates through (a) one and (b) two gluons. Graph (a) is 
impossible in the colour singlet contribution. 



of the constituents will ensure the even distribution of the hard momenta amongst the 
final constituents. The alternative is to let the constituent gluon be part of the hard 
perturbative part, that is to let it enter and end in Th- Fig. |4] (b) shows this alternative. 
Comparing Fig. |] (a) and (b), it is clear that the first is down by 1/M XJ assuming 
hard momenta run through all connecting lines in Th as is usual in SHSA. Moreover, by 
including a higher Fock state in one of the outgoing hadrons will lead to extra suppression 
due to the smaller probability of that state, not to mention the undesirable involvement 
of further unknown wavefunction. The second possibility will be our choice for colour 
neutralization in the rest of this paper. 

As we mentioned earlier, C-parity only permits the cc to annihilate through two gluons 
in the colour singlet contribution. On the contrary in the colour octet, the very presence of 
a constituent gluon in the initial state enables many more possibilities. The heavy quark- 
antiquark pair, which is now protected from the C-parity constraint by the constituent 
gluon, can annihilate through one, two and three gluons. Fig. [| (b) is an example of the 
latter and Fig. ^| (a) and (b) show a one-gluon and a two-gluon possibility, respectively. 
Fig. H] (a) and Fig. || (b) are not permitted in the singlet case. In all, the graphs can be 
divided into ten groups based on ten basic graphs from which the members of each group 
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Figure 6: The feynman graphs that represent the ten basic groups from which the graphs 
of the colour octet contributions are generated. Each graph is assigned a group number 
as indicated. 



can be constructed. These basic graphs are shown in Fig. []. The four graphs in Fig. ^| are 
in fact counted as one basic graph of group (3) since they are related by symmetry. The 
members of all groups are generated by using the afore-mentioned colour neutralization 
method. That is to let the constituent gluon end in all possible allowed position on the 
basic graph as part of T#. This generates between four to eleven graphs depending on the 
group. For example, Fig. | (b) is generated from group (5) by adding a gluon emerging 
from the \J an d ending on one of the light quark lines, Fig. [5] (a) is from group (10) 
and Fig. |5] (b) is from group (4) of Fig. || The latter two are generated by letting the 
constituent gluon to end on the three-gluon vertex on the group (10) and (4) basic graph, 
respectively. More details on these are given elsewhere |H5| . 



The decay amplitude has the form of Eq. (|3|). The convolution now involves only 
distribution amplitudes. That of the proton is shown in Sec. [|. The relevant part is 
obtained by going to transverse position 6-space and setting b = in the wavefunction 
or equivalently one can integrate out the transverse momenta k± in momentum space as 
shown in Eq. (HI). The other distribution amplitudes or wavefunctions that we need is 
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the colour octet of the xj- These have to be constructed because the importance of the 
colour octet in exclusive reactions has not been realized up to now. For xj charmonia 
with momentum P, these can be expressed as 

\X { £,P) = f /£> J [dz]45(*i, ^> *s) S$(P) \ccg; P) (27) 

where Sf^(P) is the spin-wavefunction JT^, |T^|. It is most convenient to use a product 
of delta functions for the distribution amplitude 0^ 8 ](zi, 22, 23) to set z = 0.15 for the 
constituent gluon and dividing the remaining momentum fractions Z\ = Z2 = (l — z)/2 for 
the heavy quarks. This value of z was found to be most appropriate in [17]] and the decay 
constant /® = 0.90 x 10~ 3 GeV 2 was fitted from X2 — > The same decay mode is 
not permitted for xi an d the value of is therefore unconstrained but it is expected to 
be of the same order. We found a somewhat smaller value at = 0.225 x 10~ 3 GeV 2 
for Xi from the decay into pp. 

The hard part T# of the singlet contribution is essentially the b = b' = version 
of Eq. (^) whereas that of the octet contribution now comes from an evaluation and 
summation over all members of the ten contributing groups. More explicitly the octet 
hard perturbative part can be written as 



{T^ s \x,y)} XlMM = 1 (4vra s (m c )) 3 vW o/t (2m c ) 7 



x S 9 P 9m(x,y){N^ m (x,y)} Xl ,x 2 ,x 3 

g^Groups 
m«E M ember s 



(28) 



where P gm (x,y) is the product of propagators of the member m graph of group g and 
S g is a symmetry factor of group g to take care of similar potential groups that are 
related to this group by simple change of variables. The {Ng m (x, y)}xi,x 2 M contains 
the light constituent quark helicity dependence of the outgoing nucleon and it denotes 
the numerator of this member m graph. The convolution of T# with the distribution 
amplitudes can be expressed in the form of decay form factors for the octet contributions 



B P j i8) = E 9 %^ ,4 / t da « {T^ i8 \x,y)} XlMM ||^(x,0)^(y,0)|| Al)A2)A3 , 

Ai,A 2 ,A 3 =± Amc \ qm c) J 

(29) 



where a\ = l/y/2 and o"2 = 1 as given below Eq. ([24]), and the product of proton wave- 
function factor ||\l/ p \l/ p || Al \ 2 \ 3 is similar in form to the wavefunction factor enclosed within 
braces in Eq. ( |24| ) but this one has helicity dependence and is evaluated at b = b' = 0. 



The hard part Tfj of individual feynman diagrams can be readily calculated using the 
usual feynman rules and the gauge choice by the same name. For example, the numerators 
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and propagator denominators of our example member graphs shown in Fig. |] (b), Fig. || 
(a) and (b) for some chosen helicity combinations are 



X 



6 {(z 1 -x 1 )(z 1 -y 1 )-l/A}(2m c ) 2 
1 

{(z 2 + z - x 3 )(z 2 + z-y 3 )- 1/4} (2m c ) 2 
1 1 1 



xiyi{2m c ) 2 x 2 y 2 (2m c ) 2 z(x 3 - z)(2m c 



1 



(x 3 - z)(y 3 - z)(2m c ) 2 + p 2 



(30) 



Ao.Fig. 5 (a)(^,y){iVi 7 0)Flg .| (a) (x,y)} + 



X 



(1 - xi - z)(l - yi - z){2m c ) 2 + p 2 
1 1 



(1 - z) 2 (2m c ) 2 (1 - z){\ - Vl - z)(2m c ) 



2 



x 1 11 (31) 

x 2 y 2 (2m c ) 2 x 3 y 3 (2m c ) 2 (2m c ) 2 



and 



5 (b)(^, 



5y/2 (-l) J (x 2 -y 2 ) 

■ ; "' + - + 3 (l-x 2 -*)(l-y 2 -z)(2m c ) 2 + p 2 

1 1 
{(z 2 - x 2 )(z 2 - y 2 ) - l/4}(2m c ) 2 x l2/l (2m c ) 2 

1 1 _J_ , (32 ) 

x 2 y 2 (2m c ) 2 x 3 y 3 (2m c ) 2 (2m N ' 



2 



respectively. To save on typing, we dropped +ie in the propagators but its presence should 
be implicitly understood. The p 2 in some of the denominators above does not come with 
the propagators but has been inserted by hand. The purpose of this will be explained 
below. The explicit values of S g and the form of all other P gm (x,y), {Ng m (x, 2/)}ai,a 2 ,a 3 
and ||^/ p \I/ p ||a 1 ,A2,A3 can be found in |L9 



The above expressions are merely those of three example graphs. The full calculation is 
tedious but fortunately all necessary algebra can be handled by computer. Our treatment 
and method of dealing with the colour octet contribution are largely based on and parallel 
to those in JT7|, [18] and the various parameters obtained there are used in the present 
calculations. In particular a s is set at the usual scale m c and in any case it is tied to the 
fitted octet decay constant of Xj> an d h was assumed a s ° h = ir for this special coupling 
of the constituent gluon to the perturbative part as in |]17| . This coupling is supposed 



to be softer than the others and therefore requires special treatment. Since the factors 
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' \ - v "s° ft a l wa y s appear together both in \j — > n7T and the present xj — > VP, the 
detail of each individual factor is therefore less important than the combination. If the 
calculation were done within MHSA, both a s s oit and its value need not be worried. The 
latter would be fixed dynamically. Once the constituent gluon has been dealt with and 
put into the perturbative hard part, the calculation itself is essentially the same as other 
calculations using SHSA scheme. One calculates T# perturbatively from all contributing 
graphs and then convolutes it with the distribution amplitudes of the hadron involved to 
get the probability amplitude. 

We now return to the p 2 in Eqs. (|30|), ( |3~T1) and (^) above. The purpose of inserting 
p 2 is to prevent two poles within one propagator from occurring simultaneously which is 
possible in the convolution integral in Eq. (P9f). When that happens, the ie prescription 
cannot handle such singularities. For the expression in Eq. (|30| ) the problem occurs when 
£3 = 2/3 = 2 and for Eqs. (|3l"l) and (^) this is at X\ — y\ — 1 — z and xi = 1/2 = 1 — z, 
respectively. The presence of p 2 removes this kind of singularities in the integral. The 
reason that they are there is because of the collinear approximation of SHSA. Within this 
approximation, the hadron internal transverse momentum dependence of the perturbative 
hard part is dropped T H (x, y, k±) — > T H (x, y) in favor of hard momenta running through 
the internal lines of the relevant feynman graphs. As a consequence, they are dropped 
from the denominators of the propagators as well under the assumption that the large 
virtualities of the virtual propagating particles dwarf any such dependence. Had we relin- 
quished our hold of the collinear approximation and calculated within the MHSA, such 
singularities simply do not exist. In the MHSA scheme, the internal transverse momenta 
of the hadrons involved are kept everywhere and they therefore serve as a natural cutoff 
of such singularities. From this example, the use of p 2 = (fcj_ 2 ), the average transverse 
momentum square of the nucleon, as a regulator of these singularities naturally suggested 
itself. Although this treatment is somewhat ad hoc within SHSA, it is nevertheless prac- 
tical and indeed has been used in |L7|. 

Next there is the issue of gauge invariance. The common way of expressing hadron 
wavefunctions in terms of distribution amplitudes means that the energy-momentum of 
the hadron must be divided into portions and assigned to its partonic constituents. In 
the case of charmonium in a colour octet, the energy-momentum is therefore assigned 
as P = P c + P 5 + P g with P c = ZiP, P 5 = Z2P and P g = zP. Since z\ = Z2 the mass 
of the constituent heavy quark is then Z\M XJ ~ Zi(2m c ) whereas feynman rules do not 
take into account of particles existing in a bounded, confined state and therefore give the 
standard heavy fermion propagator 1/(7 • P c — m c ) with the usual heavy quark mass m c . 
Remembering that the probability amplitude of an interaction is only gauge invariant, 
after all feynman graphs of a given order have been included, provided that the external 
particles are all on their respective mass-shells. The conflict of z\{2m c ) 7^ m c thus becomes 
a source for the violation of gauge invariance. This only arises in the graphs of group 
(3), (4) and (5). If one checks for gauge invariance, after many cancellations, the few 
remaining gauge violating terms are all proportional to (1 — 2z\)m c = zm c as expected. 
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Combining with other factors, these terms first enter at 0(z 3 ) in the numerators. Our 
calculations are therefore gauge invariant up to 0(z 2 ). As mentioned in ]nj, one could 
choose to make the results completely gauge invariant by using (1 — z)m c in the propagator 
instead of m c but because of z = 0.15, these gauge-invariance violating contributions are 
small and therefore do not worry us too much. 

Our results are shown in Table [|. The colour singlet-octet combined results are given 
along with those from the singlet alone. The singlet results here differ from those in Table 
H] because of the somewhat different values of the octet decay constants and the different 
treatments of the coupling and details between the two calculational schemes. But we are 
satisfied by the fact that the results are of the same order of magnitude. For experimental 
errors, the readers can consult ref. []3D and for theoretical ones, they will not be presented 
at this stage since they are tied closely to the experimental ones because the octet decay 
constants were obtained in |T^, [18] by the requirement of a reasonable agreement with 
two different sets of experimental measurements. As mentioned before, the disagreement 
in the measurements between PDG and BES collaboration forces us to settle around the 
average of the two sets of results. This is also the case for the decay mode xj — ¥ 71-71 '• 
Once a better agreement amongst the experiments is reached, the decay constants will 
have to be adjusted accordingly and estimates of theoretical errors be made. 

As we mentioned already in Sec. |5], in the case that the true values are close to either 
the PDG or the BES set. Colour singlet alone cannot account for the measured decay 
widths. We have outlined and described our method to include the colour octet. As can be 
seen from above, with the colour octet contribution the theoretical situation is in a much 
better shape. Note that the xj singlet wavefunction is well known and in the MHSA, 
there is no free parameter in the calculation. As to the colour octet wavefunction, that 
of the X2 is completely fixed by X2 — ► kit and therefore there is no remaining freedom 
left in y 2 — > pp. The value of /® although cannot be fixed by the pseudo-scalar decay 
mode like the X2, the smaller but still at the same order value as that of the X2 can be 
viewed as a consistent result. Since xi is a somewhat different P-wave charmonium from 
its even spin partners, one is not too surprising that the result is as it is. 

Before summarizing, we would like to address the issue of the validity of our perturba- 
tive calculations. Although the charmonium mass M XJ ~ 3.5 GeV is not extremely large, 
the scales of the coupling are set essentially by the virtualities of the intermediate gluons 
and not, for example, by the energies of the final partons. For annihilation into three 
gluons, each gluon will have about 1.1 GeV and even with a large value for Aqcd = 0.25 
GeV, the coupling will be about a s = 0.508. It is not too bad for a perturbative calcula- 
tion. On the other hand, if one compares the combined colour singlet and octet results in 
Table |] with those of the colour singlet alone in Table [|, one might be alarmed with the 
apparently large combined results when compared to the smaller singlet results and start 
questioning the validity of the perturbative calculations. First, to compare the singlet 
widths with the combined singlet and octet widths directly is misleading. A more proper 
comparison is between the singlet amplitudes with the combined amplitudes or equiva- 
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lently one compares the square root of the widths. If this is done, it can be established 
that the combined amplitudes are only a factor of 4.5 for xi an d 3.0 for \2 larger than the 
corresponding singlet ones, so the octet parts are nowhere near completely dominating 
the results. On the contrary, both the singlet and octet parts contribute comparably to 
the final results. Second, with proton-antiproton in the final states and not pions as in 
T7|j , one might worry about getting too close to the boundary of the perturbative region. 



Fortunately, the calculation of the colour singlet contribution in Sec. |5] and of J/ip decay 
into baryon-antibaryon pair in ref. J33| both within the MHSA showed that the bulk of 
the contribution came not from large coupling non-perturbative region. This should serve 
as an example of even with a pair of baryons in the final states, the calculational scheme 
of SHSA and MHSA are still under perturbative control and can be used self-consistently. 
Third, although it is not possible to do a full calculation in the MHSA, not all ten groups 
of Feynman diagrams in the colour octet contribution lead to unmanageably high dimen- 
sional integrations. For those that can be calculated numerically with sufficient accuracy, 
we have checked that the bulk of the contributions did indeed come from smaller values 
of the coupling below a s = 0.5. Using a cutoff of a c s ntoS = 0.5, well over 50 % of the 
contributions were included and with a™ toff = 0.6, 98 % were included. Of course, we 
have not confirmed reasonable widths within the MHSA due to the complexity of the 
calculations and the limitation of present days computing power, but this should still 
serve to demonstrate to a certain extent that the theoretical arguments presented above 
are reasonably sound and that the calculational schemes used are self-consistent. 

To summarize, we have provided arguments in support of the inclusion of colour 
octet in P-wave charmonium decay and by explicit calculation, have shown that although 
the same infrared divergence in inclusive decay does not exist in exclusive process, for 
consistency this next higher Fock state of the P-wave charmonium must still be included. 
With that, the would-be contradiction discussed at the end of Sec. || is resolved. The 
applicability of our arguments is obviously not restricted only to charmonia. In fact, they 
are true and should work even better for heavier quarkonia such as bb. 
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Figure Captions 



Fig. ^ The leading decay process for (a) xo an d X2 and (b) xi- A possible decay of the 
higher Fock state, the colour octet, of the xj is shown in (c). In the infrared limit 
of the outgoing gluon, (b) and (c) are degenerate. 

Fig. Feynman graphs for the colour singlet xj decay into proton-antiproton. 

Fig. ^ Percentage of the colour singlet contribution to the xj — > PP decay width that 
comes from region of values of a s below some critical cutoff value a"' 1 *'. The solid 
and dot-dashed line are for %i an d X2, respectively. This shows that the bulk of the 
contributions in fact comes from a s < 0.6 region. 

Fig. (| Example of two ways to achieve colour neutralization: (a) the constituent gluon 
enters the nucleon as part of the constituent in the nucleon after exchanging a gluon 
with the other constituents, (b) this gluon ends in the perturbative hard part T#. 

Fig. [| Example of cc annihilates through (a) one and (b) two gluons. Graph (a) is impos- 
sible in the colour singlet contribution. 

Fig. § The feynman graphs that represent the ten basic groups from which the graphs of 
the colour octet contributions are generated. Each graph is assigned a group number 
as indicated. 



Table Captions 

Table |I| The mass dimensions of the decay constants of the hadrons discussed in the text. 
They can be deduced from Eqs. (|) and (|5|). 

Table § The make-up of the parity V and C-parity C of each colour singlet and octet state 
of the J/ip and xj f° r states with up to three constituents. 

Table § Comparing colour singlet contribution of xj decay into pp to the experimental data. 

Table |] The colour singlet and octet combined partial decay width within the SHSA of 
Xj decay into proton-antiproton. Experimental measurements are also shown for 
comparison. 
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Tables 



Table 1: The mass dimensions of the decay constants of the hadrons discussed in the text. 



They can be deduced from Eqs. and (H). 



Hadron 


No. of 


Total Orbital Angular 


Mass dimension 




Constituents 


Momentum L 


of / [f] 


7T, K 


2 





1 




2 





1 


X? 


2 


1 


2 


p, p 


3 





2 




3 


1 


3 


(8) 
Xj 


3 





2 



Table 2: The make-up of the parity V and C-parity C of each colour singlet and octet 



state of the J/ip and xj f° r states with up to three constituents. 



Hadron 


LI 


Sec 


T CC 

L 9 


N 9 




C = [-l) N 9+Scc+L% 


(cc state) 














W 1J ( 3 Si) 





1 








-1 


-1 




1 


1 





1 


-1 


-1 


J/^ {8) ('So) 








1 


1 


-1 


-1 


J/ij)W ( 3 S 1 )? 





1 





1 


+1? 


+1? 


J/^ (8) ( 3 £l)? 





1 


1 


1 


-1 


+1? 


x^ ( 3 Pj) 


1 


1 








+1 


+1 







1 





1 


+1 


+1 
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Table 3: Comparing colour singlet contribution of xj decay into pp to the experimental 
data. 



J 


rW(xJ-ZP) (eV) 


PDG (eV) H 


BES (eV) @ 


1 


2.53 


75.68 ± 10.50 


37.84 


2 


16.58 


200.00 ± 20.00 


118.00 



Table 4: The colour singlet and octet combined partial decay width within the SHSA of xj 
decay into proton-antiproton. Experimental measurements are also shown for comparison. 





(eV) 


r W+W (eV) 


PDG (eV) H 


BES (eV) [0 


Xi - 


— > pp 


3.15 


56.27 


75.68 ± 10.50 


37.84 


X2 - 


— > pp 


12.29 


154.19 


200.00 ± 20.00 


118.00 
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